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Abstract 

The effect of the magnetic skew on the Parker instability is investigated by 
means of the linear stability analysis for a gravitationally stratified gas layer per- 
meated by a horizontal magnetic field. When the magnetic field is skewed (i.e., the 
field line direction is a function of the height), the wavelength of the most unstable 
mode is A ~ 10H where H is the pressure scale height. The growth rate of the 
short wavelength modes is greatly reduced when the gradient in the magnetic field 
direction exceeds 0.5 radian per scale height. Our results indicate that the Parker 
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instability in a skewed magnetic field preferentially forms large scale structures like 
giant molecular clouds. 

Subject Headings: hydromagnetics — instabilities — - interstellar: magnetic fields. 

1. INTRODUCTION 

The Parker instability has been thought to play a major role on the forma- 
tion of giant molecular clouds in galaxies (Parker 1966; Shu 1974; Mouschovias 1974; 
Mouschovias, Shu, and Woodward 1974). Although the Parker instability has many 
characteristics favorable for the formation of giant molecular clouds, it has one un- 
favorable characteristic in this regard: the fastest growing mode of the Parker in- 
stability has an infinite wave number perpendicular to the magnetic field line and 
may produce many small gas condensations instead of large scale gas condensations 
(Asseo et al. 1978). This argument has been used to suggest that the pure Parker 
instability (without self-gravity) leads to chaotic structures and turbulence rather 
than large scale cloud complexes (Asseo et al. 1978; Elmegreen 1982). This ques- 
tion, however, is based on a linear stability analysis which assumes that the magnetic 
field is parallel in the unperturbed state. The Galactic magnetic field, however, is 
not perfectly parallel (Heiles 1987; Sawa & Fujimoto 1986; Shibata & Matsumoto 
1991), but rather disturbed and skewed. When the magnetic field is skewed, the 
mode having a large wave number perpendicular to the magnetic field is likely to be 
suppressed, since such modes force the adjacent field lines to intersect. Thus, the 
fastest growing mode is expected to have a finite wave number and to produce large 
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scale condensations. This Letter studies the effect of magnetic skew on the Parker 
instability by means of linear stability analysis and supports the above expectation. 



The stability analysis is carried out in §2 and its application to giant molecular 
cloud formation is discussed in §3. 



We consider an idealized model of a skewed magnetic field to evaluate the effect 
of the magnetic skew on the Parker instability. The model consists of an isothermal, 
stratified gas layer permeated by a horizontal magnetic field whose direction changes 
as a function of height at a constant pitch, <p (z) = tan -1 (B y /B x ) oc z. The 
density distribution and the magnetic field of the model are 



where H is the scale height and e denotes the strength of the magnetic skew. (The 
magnetic skew is not necessary to be constant to suppress the growth of short wave 
length modes. We imposed this assumption in order to minimize the number of model 
parameters and to make our model as simple as possible.) This model atmosphere 
is in hydrostatic equilibrium under gravity in the ^-direction. From the equation of 
motion, the scale height is expressed as 



2. LINEAR STABILITY ANALYSIS 



p = p exp (-z/H) , 



(1) 



B = (B x , B yi B z ) 



(2) 



[B cos (ez/H) exp (-z/2H), B sin (ez/H) exp (-z/2H), 0] 



gH = c s 2 + 




(3) 
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where g is the gravitational acceleration and assumed to be constant for simplicity. 
The plasma beta, (3 = 8irpc s 2 /B 2 , is constant in this model atmosphere, and is 
assumed to be unity, f3 = 1, in the following except when otherwise is noted. Figure 
1 shows a schematic view of the model magnetic field. (Although this configuration 
of the magnetic field is called "magnetic shear" in the solar physics, we use the word 
"skew" in order to avoid the confusion with the magnetic field sheared with the 
Galactic differential rotation.) 

We consider a small perturbation superimposed on the equilibrium described 
above. The perturbation is assumed to be isothermal and to be described by the 
ideal MHD equations. Since the perturbation equations are standard, their detailed 
form is omitted here (see, e.g., Horiuchi et al. 1988 for perturbation equations). A 
normal mode solution of the perturbation equations has the form of 

Pi(x, y, z) = pi_{z) exp {-iut + ik x x + ik y y) . (4) 

According to the variational principle (see, e.g., Bernstein 1983; Nakamura, Hanawa, 
and Nakano 1991) the square of the eigenfrequency, a; 2 , is real. The boundary con- 
dition is set so that the energy density of the eigenfunction is periodic, 

Pl (z + 2nH/e) = e ie ~^ e Pl (z) . (5) 

See Shu (1974) and Nakamura, Hanawa, and Nakano (1991) for the reason why the 
energy density should be constant in an eigenfunction. The ratio of the phase shift, 
6>, to the vertical period of the magnetic skew, 2nH/e, is an effective wavenumber 
in the vertical direction (/c z , e // = 9e/2irH). Since our model has the helical 

4 



symmetry, the eigenfrequency is a function of the total horizontal wave number, u> = 
^{\/k x 2 + ky 2 ). [ Our model is invariant under the transformation of {x, y, z) — > 
(xcos5 + ysin5, — xsin5 + ycos5, 2; + SH/s), where 5 is an arbitrary number. ] 

Figure 2 shows the stability diagram of our model atmosphere for e = 1 and 
/? = 1. The ordinate and the abscissa are the growth rate, —u 2 H/g : and the wave 
number, k 2 H 2 , respectively. The growth rate and the wave number are normalized 
by \fgjH and 1/H , respectively. The stability diagram has a band structure similar 
to the energy level diagram of a crystal (cf. Kittel 1976). The darkly painted regions 
denote the allowed bands and the blank regions show the inhibited bands. The most 
unstable mode has a finite wave number of kH = 0.70. 

On the top of the uppermost band v\ z of the eigenmode has no node in the 
vertical direction {9 = 0). The value of 9 changes from 2mr to 2(n + 1)tt in each band. 
As the number of nodes in the eigenfunction increases (9 increases), the growth rate 
decreases continuously but with some gaps. This is a general property of the Parker 
instability. When the initial magnetic field is parallel, the growth rate decreases 
continuously as the wavenumber, \k z \, increases (see, e.g., Parker 1979). When the 
Parker unstable layer has a finite thickness, the growth rate is discrete and decreases 
with the increase in the number of nodes in the eigenfunction (Horiuchi et al. 1988). 
In our model, the Parker unstable layers appear periodically in the vertical direction 
for a given k. In the band structure of figure 2 the growth rate is continuous in 
each band and discrete between bands. This is because our model atmosphere has a 
periodical structure where each Parker unstable layer has a finite thickness but the 
number of the Parker unstable layers is infinite. 
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Figure 3 shows the maximum growth at a given wave number for £ = 0.2, 
0.5, and 1.0. The growth rate decreases as e increases for a given k. The decrease 
is larger for a larger wave number. When e > 0.2, the mode of k = oo is not the 
fastest growing mode. The most unstable mode has a wave length of A ~ 10H for 
0.2 < e < 1.0. 

Figure 4 shows the maximum growth at k = oo as a function of s. As e 
increases, the growth rate at k = oo decreases. The growth of the mode having k 
is suppressed by magnetic tension when k ■ B is large. Whenever the magnetic field 
is not parallel and k is sufficiently large, most magnetic field lines are bent with a 
short wavelength and the magnetic tension suppresses the growth of the mode. From 
figures 3 and 4 we find that the magnetic skew suppresses the Parker instability of 
short wavelength modes and that the effect is appreciable for e > 0.5. Note that 
e = 0.5 corresponds to the case when the direction of the magnetic field changes 
29° per scale height. 

3. APPLICATION TO GIANT MOLECULAR CLOUD FORMATION 

The analysis of the previous section demonstrated that the fastest growing 
Parker instability has a wavelength of A = 10H when the magnetic field is skewed 
appreciably. It implies that large scale mass condensations are formed by the Parker 
instability if the Galactic magnetic field is skewed significantly, i.e., more than several 
tens degrees per scale height. The Galactic magnetic field is not uniform in the Solar 
neighborhood (see, e.g., Heiles 1987 and the references therein) and has a substantial 
nonuniform component. Galactic dynamo theory (see, e.g., Sawa and Fujimoto 1986; 
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Fujimoto and Sawa 1987) also predicts that the Galactic magnetic field changes its 
direction as a function of the height; in their model the direction of the magnetic 
field in the halo (at the height of z = 1 kpc ~ 6H) differs 180° from that 
in the disk (e ~ 0.5). Also the nonlinear development of the Parker instability 
produces the skewed magnetic field owing to the Coriolis force, even when the initial 
field is not skewed (Shibata and Matsumoto 1991). Consequently, although a firm 
observational evidence of the magnetic skew has not been obtained, it is very likely 
that the Galactic magnetic field has a skew component which is significantly strong 
for the formation of giant molecular clouds through the Parker instability. 

Our results (Fig. 3) indicate that the growth time of the Parker instability is 

7 f u 2 H/g \-i/2 / H \l/2 / g x-i/2 

r ~ 4.9xl0 7 ^ =■) years, 6) 

V5.8xl0" 2 / V160pc/ V3.5 x 10" 9 cms" 2 / y KI 

in a skewed field with e ~ 0.2 — 1.0 and f3 = P g /P m = 1 and the mass of the 
cloud formed by the Parker instability is 

M ~ 2p\ 2 H ~ 2 x 10 7 ( ( —Y ( — ^— ) 3 M . (7) 
1 VI cm" 3 / \10HJ Vl60pc/ ° V ; 

We now briefly comment on the effect of the self-gravity on the Parker insta- 
bility. Hanawa, Nakamura, and Nakano (1991) has analyzed the linear stability of 
the Galactic gaseous disk taking account of the magnetic field, the Galactic rotation, 
and the self-gravity of the gas. The magnetic field was assumed to be parallel in 
the equilibrium. According to them, the growth rate is larger when the wavenumber 
perpendicular to the field line is larger. The growth of short wavelength modes are 
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not suppressed by the self-gravity. "Turbulence" rather than a cloud is likely to be 
formed as far as the magnetic field is parallel in equilibrium even if the self-gravity 
is taken into account. Thus, the skewed magnetic field is an important factor to 
produce a large scale structure like clouds through the Parker instability. 

The skewed magnetic field may affect not only the linear growth rate of the 
Parker instability but also the nonlinear evolution of the Parker unstable modes. 
The intersection of the magnetic field lines is strictly inhibited in the framework 
of the ideal magnetohydrodynamics. Thus, the growth of an unstable mode will 
be saturated at latest at the stage when simple extrapolation of the linear stability 
analysis predicts the intersection of the magnetic field lines (see Matsumoto et al 
1988, 1990 for justification of this conjecture). Figure 5 shows the eigenfunction of 
the Parker unstable mode in the skewed magnetic field (e = 0.5 and ,9 = 1). The 
ordinate is the velocity perturbation, v± z , at the origin ( x, y) = (0, 0) and the 
abscissa is the height z. The thick curve is for the mode of (/c x , k y ) = ( 0.7/if, 0) 
and the thin curve for that of (k x , k y ) = ( 0, 0.7 /H). As shown in the figure the 
velocity, vi jZ , is not a monotonically increasing function of z. At some height, the 
rising magnetic field line catches up with the foregoing field line as the perturbation 
grows. The growth of the perturbation is likely to be saturated as the field lines 
becomes closer to each other. If the two eigenmodes shown in figure 5 are excited 
simultaneously, however, the field lines are less likely to catch up with the neighboring 
field lines. The linear combination of the eigenmodes is also an eigenmode since the 
growth rates are the same [w(k) = cu(|k|)]. It implies that the perturbation is 
less easy to be saturated when more than 2 plane waves are simultaneously excited. 
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If all the mode having the same wave number, |k|, are excited, the perturbation 
velocity increases with height monotonically [v\ z oc Pq~ 1 ^ 2 oc exp (z/2H)] at 
( x, y) = ( 0, 0). Then, the perturbation grows up to have a large amplitude 
without suffering serious nonlinear effects. The linear combination of many plane 
waves leads to a cylindrical wave having a symmetry axis at ( x, y) = ( 0, 0) 
(see Fig. 5). Thus, we expect that a circular structure would be formed in a skewed 
magnetic field rather than a stripe. 
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Figure Legends 

Fig. 1. Schematic view of the skewed magnetic field. 

Fig. 2. The stability diagram for e = 0.5 and (3=1. The ordinate is the square of 
the growth rate, — (H/g)uj 2 , and the abscissa is the square of the wave number, 
H 2 k 2 . The eigenmode has a band structure. The allowed bands are denoted 
with black and the inhibited bands are denoted with blank. 

Fig. 3. The maximum growth rate at a given wave number for e = 0.2, 0.5, and 1.0. 

Fig. 4. The maximum growth rate at k = oo as a function of e. 

Fig. 5. The velocity perturbation, v\ z , of the unstable eigenmode at ( x, y) = ( 0, 0) 
as a function of the height, z for s = 0.5, f3 = 1.0, |k|if = 0.7, and uj 2 H/g = 
—5.82 x 10~ 2 . The thick solid curve is for the mode of ( k x H, k y H) = 
( 0.7, 0.0) and the thin solid curve for that of ( k x H, k y H) = ( 0.0, 0.7). The 
dashed curve denotes the linear combination of many plane wave eigenmodes, 
where the velocity perturbation is proportional to v\ z oc p~ x l 2 at ( x, y) = 
(0, 0). 
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